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The critical current needed to destabilize the magnetization of a perpendicular ferromagnet via 
the spin Hall effect is studied. Both the dampinglike and fieldlike torques associated with the spin 
current generated by the spin Hall effect is included in the Landau-Lifshitz-Gilbert equation to 
model the system. In the absence of the fieldlike torque, the critical current is independent of the 
damping constant and is much larger than that of conventional spin torque switching of collinear 
magnetic systems, as in magnetic tunnel junctions. With the fieldlike torque included, we find that 
the critical current scales with the damping constant as a 0 (i.e., damping independent), a , and 
<a 1//2 depending on the sign of the fieldlike torque and other parameters such as the external field. 
Numerical and analytical results show that the critical current can be significantly reduced when 
the fieldlike torque possesses the appropriate sign, i.e. when the effective field associated with the 
fieldlike torque is pointing opposite to the spin direction of the incoming electrons. These results 
provide a pathway to reducing the current needed to switch magnetization using the spin Hall effect. 

PACS numbers: 75.78.-n, 75.70.Tj, 75.76.+j, 75.40.Mg 


I. INTRODUCTION 

The spin Hall effect^ - — (SHE) in a nonmagnetic heavy 
metal generates pure spin current flowing along the di¬ 
rection perpendicular to an electric current. The spin 
current excites magnetization dynamics in a ferromagnet 
attached to the nonmagnetic heavy metal by the spin- 
transfer effect^. There have been a number of exper¬ 
imental reports on magnetization switching and steady 
precession induced by the spin Hall effect^ - —. These dy¬ 
namics have attracted great attention recently from the 
viewpoints of both fundamental physics and practical ap¬ 
plications. 

An important issue to be solved on the magnetization 
dynamics triggered by the spin Hall effect is the reduc¬ 
tion of the critical current density needed to destabilize 
the magnetization from its equilibrium direction, which 
determines the current needed to switch the magneti¬ 
zation direction or to induce magnetization oscillation. 
The reported critical current density for switching^iiS - — 
or precession^ is relatively high, typically larger than 10 7 
A/cm 2 . One of the reasons behind this may be related 
to the recently predicted damping constant independent 
critical current when SHE is use d 14 ! 15 . This is in con¬ 
trast to spin-transfer-induced magnetization switching in 
a typical giant magnetoresistance (GMR) or magnetic 
tunnel junction (MTJ) device where the critical current 
is expected to be proportional to the Gilbert damping 
constant a. Here the magnetization dynamics is excited 
as a result of the competition between the spin torque 
and the damping torquei^. Since the damping constant 
for typical ferromagnet in GMR or MTJ devices is rela¬ 
tively small (a ~ 10 -2 — 10~ 3 ) 17 d 8 , it can explain why 
the critical current is larger for the SHE driven systems. 
Thus in particular for device application purposes, it is 
crucial to find experimental conditions in which the mag¬ 


netization dynamics can be excited with lower current. 

Another factor that might contribute to the reduc¬ 
tion of the critical current is the presence of the field 
like torque^. In the GMR/MTJ systems, both the con¬ 
ventional spin torque, often referred to as the damp¬ 
inglike torque, and the fieldlike torque arise from the 
spin transfer between the conduction electrons and the 
magnetizat iorAi^r— . Due to the short relaxation length 
of the transverse spin of the conduction electron s 24 ! 25 , 
the damping like torque is typically larger than the field¬ 
like torque. Indeed, the magnitude of the field like 
torque experimentally found in GMR/MTJ systems has 
been reported to be much smaller than the damping like 
torque^""—. Because of its smallness, the fieldlike torque 
had not been considered in estimating the critical current 
in the GMR/MTJ systemsi&22r— , although it does play a 
key role in particular system s 33 ! 34 . In contrast, recent ex¬ 
periments found that the fieldlike torque associated with 
the SHE is larger than the damping like torque^ - — . 

The physical origin of the large SHE-induced field¬ 
like torque still remains unclear. Other possible sources 
can be the Rashba effec t 36 1 41 ~ —, bulk effect^, and the 
out of plane spin orbit torque^. Interestingly, the field 
like torque has been reported to show a large angu¬ 
lar dependenc e 36 ! 37 ! 47 (the angle between the current 
and the magnetization), which cannot be explained by 
the conventional formalism of spin-transfer torque in 
GMR/MTJ systems. The fieldlike torque acts as a torque 
due to an external field and modifies the energy landscape 
of the magnetization. As a result, a large fieldlike torque 
can significantly influence the critical current. However, 
the fieldlike torque had not been taken into account in 
considering the current needed to destabilize the magne¬ 
tization from its equilibrium direction and thus its role 
is still unclear. 

In this paper, we study the critical current needed to 
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destabilize a perpendicular ferromagnet by the spin Hall 
effect. The Landau-Lifshitz-Gilbert (LLG) equation with 
the dampinglike and fieldlike torques associated with the 
spin Hall effect is solved both numerically and analyti¬ 
cally. We find that the critical current can be significantly 
reduced when the fieldlike torque possesses the appropri¬ 
ate sign with respect to the dampinglike torque. With 
the fieldlike torque included, the critical current scales 
with the damping constant as <a° (i.e., damping indepen¬ 
dent), a, and ad/ 2 , depending on the sign of the field- 
like torque and other parameters. Analytical formulas 
of such damping-dependent critical current are derived 
[Eqs. (USD-®], and they show good agreement with the 
numerical calculations. From these results, we find con¬ 
ditions in which the critical current can be significantly 
reduced compared to the damping-independent thresh¬ 
old, i.e., systems without the fieldlike torque. 

The paper is organized as follows. In Sec. im we 
schematically describe the system under consideration. 
We discuss the definition of the critical current in Sec. 
uni Section 113 summarizes the dependences of the crit¬ 
ical current on the direction of the damping constant, 
the in-plane field, and the fieldlike torque obtained by 
the numerical simulation. The analytical formulas of the 
critical current and their comparison to the numerical 
simulations are discussed in Sec. El The condition at 
which damping-dependent critical current occurs is also 
discussed in this section. The conclusion follows in Sec. 

ED 


II. SYSTEM DESCRIPTION 

The system we consider is schematically shown in Fig. 
[lj where an electric current flowing along the x-direction 
injects a spin current into the ferromagnet by the spin 
Hall effect. The magnetization dynamics in the ferro¬ 
magnet is described by the LLG equation, 

dm dm 

dt 1 dt (1) 

— jH s m x (e y x m) — 7/3iL s m x e y , 

where 7 and a are the gyromagnetic ratio and the Gilbert 
damping constant, respectively. We assume that the 
magnetization of the ferromagnet points along the film 
normal (i.e., along the z axis), and an external in-plane 
magnetic field is applied along the x or y axis. The total 
magnetic field H is given by 

H = i?appi n H + H K m z e z , ( 2 ) 

where H app i is the external field directed along the x or 
y axis and Hk is the uniaxial anisotropy field along the 
z axis, n# and e* are unit vectors that dictate the di¬ 
rection of the uniaxial anisotropy field and the i axis, 
respectively. Here we call the external field along the x 
and y directions the longitudinal and transverse fields, 
respectively. The third and fourth terms on the right- 
hand side of Eq. (pQ) are the damping like and fieldlike 



FIG. 1. Schematic view of the spin-Hall system. The x 
axis is parallel to current, whereas the z axis is normal to the 
film plane. The spin direction of the electrons entering the 
magnetic layer via the spin Hall effect points along the +y or 
—y direction. 


torques associated with the spin Hall effect, respectively. 
The torque strength H s can be expressed with the current 
density j, the spin Hall angle $, the saturation magneti¬ 
zation M, and the thickness of the ferromagnet d, i.e., 


2 eMd 


( 3 ) 


The ratio of the fieldlike torque to the damping like 
torque is represented by / 3 . Recent experiments found 
that f 3 is positive and is larger than 1— . 

The magnetization dynamics described by the LLG 
equation can be regarded as a motion of a point particle 
on a two-dimensional energy landscape. In the presence 
of the fieldlike torque, the energy map is determined by 
the energy density given by^i 


£ = —M 


/ 


dm • H — / 3 MH s m • e y . 


( 4 ) 


Then, the external field torque and the fieldlike torque, 
which are the first and fourth terms on the right-hand- 
side of Eq. ©, can be expressed as —7m x B , where the 
effective field B is 


d£ 

d Mm 


( 5 ) 


The initial state of the numerical simulation is chosen to 
be the direction corresponding to the minimum of the 
effective energy density £. The explicit forms of the ini¬ 
tial state for the longitudinal and the transverse external 
fields are shown in Appendix A. 

We emphasize for the latter discussion in Sec. IVIthat. 
using Eqs. ©, ®, and ©, the time change of the effec¬ 
tive energy density is described as 


d£ d£ s d£ a 
dt dt ^ dt 


( 6 ) 








3 


Here the first and second terms on the right-hand side 
are the rates of the work done by the spin Hall torque 
and the dissipation due to damping, respectively, which 
are explicitly given by 

-jjjT — lMH s [ej, • B — (m • e y ) (m • B )\, (7) 


d£ a 

dt 


—a^M 


B 2 


(m • B ) 2 


(8) 


The sign of Eq. © depends on the current direction 
and the effective magnetic field, while that of Eq. m is 
always negative. 

The magnetic parameters used in this paper mimic the 
conditions achieved in CoFeB/MgO heterostructures^; 
M = 1500 emu/c.c., Hk = 540 Oe, d = 0.1, 7 = 
1.76 x 10 7 rad/(Oe s), and d = 1.0 nm. The value of 
/3 is varied from — 2 , 0, to 2 . Note that we have used a 
reduced (Refs^^) in order to obtain critical currents 
that are the same order of magnitude with that obtained 
experimentally. We confirmed that the following discus¬ 
sions are applicable for a large value of Hk IT). 


III. DEFINITION OF CRITICAL CURRENT 

In this section, we describe how we determine the crit¬ 
ical current from the numerical simulations. In exper¬ 
iments, the critical current is determined from the ob¬ 
servation of the magnetization reversa l 12 i 41 ! 46 i 48 ~ —. As 
mentioned in Sec. nn in this paper, the initial state for 
calculation is chosen to be the minimum of the effective 
energy density. Usually, there are two minimum points 
above and below the xy plane because of the symmetry. 
Throughout this paper, the initial state is chosen to be 
the minimum point above the xy plane, i.e., m z ( 0 ) > 0 , 
for convention.” It should be noted that, once the mag¬ 
netization arrives at the xy plane during the current ap¬ 
plication, it can move to the other hemisphere after the 
current is turned off due to, for example, thermal fluc¬ 
tuation. Therefore, here we define the critical current as 
the minimum current satisfying the condition 

lim m z (t) < e, (9) 

t—yoo 

where a small positive real number e is chosen to be 
0.001. The duration of the simulations is fixed to 5 /is, 
long enough such that all the transient effects due to the 
current application are relaxed. Figures [ 2 ja) and [ 2 jb) 
show examples of the magnetization dynamics close to 
the critical current, which are obtained from the numer¬ 
ical simulation of Eq. ©• As shown, the magnetization 
stays near the initial state for j = 3.1 x 10 6 A/cm 2 , while 
it moves to the xy plane for j = 3.2 x 10 6 A/cm 2 . Thus, 
the critical current is determined as 3.2 x 10 6 A/cm 2 in 
this case. 

We note that the choice of the definition of the criti¬ 
cal current has some arbitrariness. For comparison, we 
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FIG. 2. Time evolution of the z component of the magnetiza¬ 
tion m z in the presence of the transverse field of H app \ = 200 
with (a) j = 3.1 x 10 6 A/cm 2 and (b) j = 3.2 x 10 6 A/cm 2 . 
The value of f5 is zero. 


show numerically evaluated critical current with a differ¬ 
ent definition in Appendix B. The main results of this 
paper, e.g., the dependence of the critical current on the 
damping constant, are not affected by the definition. 

We also point out that the critical current defined by 
Eq. (j9]) focuses on the instability threshold, and does 
not guarantee a deterministic reversal. For example, 
in the case of Fig. [2fb), the reversal becomes prob¬ 
abilistic because the magnetization, starting along +z, 
stops its dynamics at the xy plane and can move back 
to its original direction or rotate to a point along —2: 
resulting in magnetization reversal. Such probabilistic 
reversal can be measured experimentally using transport 
measurements ^ 12 ! 41 ! 46 ! 49 ! 50 or by studying nucleation of 
magnetic domains via magnetic imaging^. On the other 
hand, it has been reported that deterministic reversal can 
take place when a longitudinal in-plane field is applied 
alongside the curren t 41 ! 49 . It is difficult to determine the 
critical current analytically for the deterministic switch¬ 
ing for all conditions since, as in the case of Fig. [2jb), 
the magnetization often stops at the xy plane during the 
current application. This occurs especially in the pres¬ 
ence of the transverse magnetic field because all torques 
become zero at m = ±e y and the dynamics stops. Here 
we thus focus on the probabilistic reversal. 
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FIG. 3. Numerically evaluated m z at t — 5 /is for (a)-(c) the longitudinal ( 11 # = e x ) and (d)-(f) the transverse (n h = e y ) 
fields, where the value of (3 is (a), (d) 0.0; (b), (e) 2.0; and (c), (f) —2.0. The damping constant is a — 0.005. The color scale 
indicates the z component of the magnetization (m z ) at t = 5 /is. The red/white boundary indicates the critical current for 
probabilistic switching, whereas the red/blue boundary gives the critical current for deterministic switching. 


IV. NUMERICALLY ESTIMATED CRITICAL 
CURRENT 


In this section, we show numerically evaluated critical 
current for different conditions. We solve Eq. m and 
apply Eq. m to determine the critical current. Figure 
[3] shows the value of m z at t = 5 /is in the presence of 
(a)-(c) the longitudinal (n h = e^) and (d)-(f) the trans¬ 
verse (n# = e y ) fields. The value of f3 is 0 for Figs. 
(3](a) and (3](d), 2.0 for Figs. [3Kb) and [31(e), and —2.0 for 
Figs. He) and [3](f), respectively. The damping constant 
is a = 0.005. The red/white boundary indicates the crit¬ 
ical current for the probabilistic switching, whereas the 
red and blue (m z = — 1 ) boundary gives the critical cur¬ 
rent for the deterministic switching. Using these results 
and the definition of the critical current given by Eq. ©, 
and performing similar calculations for different values of 
<a, we summarize the dependence of the critical current on 
the longitudinal and transverse magnetic fields in Fig. [4j 
The damping constant is varied as the following in each 
plot: a = 0.005, 0.01, and 0.02. The solid lines in Fig. [4] 
represent the analytical formula derived in Sec. [VI 


A. In the presence of longitudinal field 

In the case of the longitudinal field and /3 = 0 
shown in Fig. 0(a), the critical current is damping- 
independent. Such damping-independent critical current 
has been reported previously for deterministic magneti¬ 
zation switchin g 14 ! 15 . Similarly, in the case of the longi¬ 
tudinal field and negative /3 (/3 = —2.0) shown in Fig. 
He), the critical current is damping-independent. In 
these cases, the magnitude of the critical current is rel¬ 
atively high. In particular, near zero field, the critical 
current exceeds ~ 10 8 A/cm 2 , which is close to the limit 
of experimentally accessible value. These results indicate 
that the use of the longitudinal field with zero or negative 
/3 is ineffective for the reduction of the critical current. 

On the other hand, when /3 is positive, the critical cur¬ 
rent depends on the damping constant, as shown in Fig. 
Hb). Note that positive /3 is reported for the torques 
associated with the spin Hall effect or Rashba effect in 
the heterostructures studied experimentally^ - —^. The 
magnitude of the critical current, ^ 10 x 10 6 A/cm 2 , is 
relatively small compared with the cases of zero or neg¬ 
ative /?. In this case, the use of a low damping material 
is effective to reduce the critical current. Interestingly, 
the critical current is not proportional to the damping 
constant, while that previously calculated for a GMR or 
MTJ system^ is proportional to a. For example, the 
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FIG. 4. Numerically evaluated critical currents in the presence of (a)-(c) the longitudinal (njr = e : , ; ) and (d)-(f) the transverse 
(n h = e, ( ) flelds, where the value of 8 is (a), (d) 0.0; (b), (e) 2.0; and (c), (f) —2.0, respectively. The solid lines are analytically 
estimated critical current in Sec. El 


critical current at zero longitudinal field in Fig. ID^b) is 
12.3, 17.2, and 24.0 x 10 6 A/cm 2 for a = 0.005, 0.01, and 
0.02, respectively. These values indicate that the critical 
current is proportional to a 1 ' 2 . In fact, the analytical 
formula derived in Sec. IVl shows that the critical current 
is proportional to a 1 / 2 for positive /3 [see Eq. (HU)]. 

To summarize the case of the longitudinal field, the 
use of a heterostructure with positive /3, which is found 
experimentally, has the possibility to reduce the critical 
current if a ferromagnet with low damping constant is 
used. In this case, the critical current is proportional to 
a 1 / 2 , which has not been found in previous works. 


B. In the presence of transverse field 

In the presence of the transverse field with /3 = 0, 
the critical current shows a complex dependence on the 
damping constant a, as shown in Fig. |U[d). When the 
current and the transverse field are both positive (or neg¬ 
ative) , the critical current is proportional to the damping 
constant a except near zero field. The numerically cal¬ 
culated critical current matches well with the analytical 
result, Eq. shown by the solid lines. In this case, 

the use of the low damping material results in the reduc¬ 
tion of the critical current. On the other hand, when the 
current and the transverse field possess the opposite sign, 
the critical current is damping independent. Moreover, in 
this case, the critical current is of the order of 10 8 A/cm 2 . 
Thus, it is preferable to use the current and field having 
the same sign for the reduction of the critical current. It 


should be noted that, in our definition, the same sign of 
current and field corresponds to the case when the direc¬ 
tion of incoming electrons’ spin (due to the SHE) and the 
transverse field are opposite to each other. The reason 
why the critical current becomes damping dependent in 
this situation will be explained in Sec. IVl 

When /3 is positive the critical current depends on the 
damping constant for the whole range of the transverse 
field, as shown in Fig. IU[e). The critical current is 
roughly proportional to ad/ 2 , in particular, close to zero 
field. The solid lines display the analytical formula, Eq. 
dm and show good agreement with the numerical calcu¬ 
lations. The damping dependence of the critical current 
becomes complex when the magnitude of the transverse 
field is increased [see Eq. dm We note that the critical 
current for the positive f3 in Fig. |U[e) is smaller than that 
for /3 = 0 in Fig. |U[d) for the whole range of H a pp i. 

On the other hand, when /3 is negative, the critical 
current is almost independent of a, especially near zero 
field. However, when the transverse field is increased, 
there is a regime where the critical current depends on 
the damping constant. Such transition of the critical 
current with the transverse field is also predicted by the 
analytical solution, Eq. <m 

To summarize the case of the transverse field, the a de¬ 
pendence of the critical current can be categorized into 
the following: a 0 (damping independent), a, a 1 / 2 , or 
other complex behavior. As with the case of the longi¬ 
tudinal field, the use of a heterostructure with positive /3 
allows reduction of the critical current when low damping 
ferromagnet is used. Overall, the most efficient condition 
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to reduce the critical current is to use the transverse field 
with heterostructures that possess low a and positive (3. 
In this case, the critical current is reduced to the order 
of 10 6 A/cm 2 . 

V. ANALYTICAL FORMULA OF CRITICAL 
CURRENT 

In this section, we derive the analytical formula of the 
critical current from the linearized LLG equation^. The 
complex dependences of the critical current on the damp¬ 
ing constant a discussed in Sec. [TV] are well explained by 
the analytical formula. We also discuss the physical in¬ 
sight obtained from the analytical formulas. 


By = i^appi sin# cos(p— ph)+/3H s sinO sin p+H^ cos 2 0, 

(17) 

where pn represents the direction of the external field 
within the xy plane: pn = 0 for the longitudinal field 
and 7r/2 for the transverse field. 

The solution of Eq. m is mx,rriY oc 
exp{ 7 [±iy // det[M] — (Tr[M]/2) 2 — Tr [M] / 2] t}, where 
det[M] and Tr[M] are the determinant and trace of 
the matrix M, respectively. The imaginary part of the 
exponent determines the oscillation frequency around 
the Z axis, whereas the real part determines the time 
evolution of the oscillation amplitude. The critical 
current is defined as the current at which the real part 
of the exponent is zero. Then, the condition Tr[M] = 0 
gives 


A. Derivation of the critical current 

To derive the critical current, we consider the stable 
condition of the magnetization near its equilibrium. It is 
convenient to introduce a new coordinate XYZ in which 
the Z axis is parallel to the equilibrium direction. The 
rotation from the xyz- coordinate to the XYZ coordinate 
is performed by the rotation matrix 

( cos# 0 — sin#\ / cos(/? sin p 0\ 

0 1 0 J | — sinp cos p 0 ] , (10) 

sin 6 0 cos 6 / \ 0 0 1 J 


a {Bx + By) — 2 H s sin 0 sin p = 0, (18) 

For the longitudinal field, Eq. m gives 

•long 2eYaMd \JpXz -^appi 

m PWTW) ’ 

indicating that the critical current is roughly propor¬ 
tional to ad/ 2 . This formula works for positive (3 only^ 
if we assume 0 < 2 + a(3 ~ 2, which is satisfied for typical 
ferromagnets. The critical current when the transverse 
field is applied reads 


where (0, p) are the polar and azimuth angles of the 
magnetization at equilibrium. The equilibrium magne¬ 
tization direction under the longitudinal and transverse 
magnetic field is given by Eqs. m and (jA2|h respec¬ 
tively. Since we are interested in small excitation of the 
magnetization around its equilibrium, we assume that 
the components of the magnetization in the XYZ coor¬ 
dinate satisfy mz — 1 and |mx|, \tuy\ <C 1. Then, the 
LLG equation is linearized as 



•TRANS 

Jc 


2 aeMd 

Wh—Jh k) k 


1 / -f^appl 

2 v~^r 


2 ' 


( 20 ) 


when (3 = 0, indicating that the critical current is pro¬ 
portional to a. The critical current for finite (3 is 


•trans eMd 

Jc M 

— (1 + Oi(3)H a p p i =b yjH 2 ppl + 2af3{2 + ot(3)H-^ 
X PV + clP) 

( 21 ) 


where the components of the 2x2 matrix M are 


Mi,i “ aBx ~ H s sin ^ sin 99, (12) 

Mi, 2 = B y , (13) 

M2,i = Bx (14) 

M2 5 2 = glBy — H s sin 0 sin p. (15) 

Here, Bx and By are defined as 


Equation m works for the whole range of |iJ a ppi|(< 
Hk) for positive /3, while it only works when |if a ppi| > 
2a/3(2 + a/3)HK for negative (3. For example, when 
/3 = —2.0, this condition is satisfied when |iJ app i| >108 
Oe for a = 0.005 and |# app i| > 152 Oe for a = 0.01. 
However the condition is not satisfied for the present 
range of H app i for a = 0.02. The solid lines in Fig. |4jf) 
show where Equation m is applicable. The zero-field 
limits of Eqs. m and m become identical, 


lim j c = ± 

^appl ^"0 


2 ey/aMd 
M 


V%Hk 
(3(2 + a/3) 


( 22 ) 


Bx = -H app i sinO cos(p—ph)-\-/3H s sin^ sin p+H^ cos 20, 

( 16 ) 


indicating that the critical current near zero field is pro¬ 
portional to ai 1 / 2 when (3 > 0. 
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time (ns) 

FIG. 5. Magnetization dynamics under the conditions of (a) 
n h — e y , H app \ — 50 Oe, = 0, a = 0.005, and j — 13.2 x 10 6 
A/cm 2 , and (b) uh — e x , H& PP i = 50 Oe, — 0, a = 0.005, 
and j — 90 x 10 6 A/cm 2 . 

B. Discussions 

The solid lines in Fig. @|b), SJd), ®[e), and|U^f) show 
the analytical formulas, Eqs. (H9]b ([20]) . and (]2T|) . As 
evident, these formulas agree well with the numerical re¬ 
sults in the regions where the critical currents depend on 
the damping constant. In this section, we discuss the rea¬ 
son why the critical current becomes damping dependent 
or damping independent depending on the field direction 
and the sign of /?. 

It is useful for the following discussion to first study 
typical magnetization dynamics found in the numerical 
calculations. Figure [5] shows the time evolution of the 
x, y and 2 components of the magnetization when the 
critical current depends on [Fig. Ef a )] or is independent 
of [Fig. [5][b)] the damping constant. For the former, 
the instability is accompanied with a precession of the 
magnetization. On the other hand, the latter shows that 
the instability takes place without the precession. 

We start with the case when the critical current be¬ 
comes damping dependent. To provide an intuitive pic¬ 
ture, we schematically show in Fig. E*a) the torques ex¬ 
erted on the magnetization during one precession period 
when current is applied. The condition is the same with 
that described in Fig. O a )> be., th e transverse magnetic 
field is applied with (3 = 0. In Fig. E*a), magnetization 
is shown by the large black arrow, while the directions 
of the spin Hall torque, the damping torque and the ex¬ 


ternal field torque are represented by the solid, dotted 
and dashed lines, respectively (the external field torque 
is tangent to the precession trajectory). As evident in 
Fig. [5][a), the precession trajectory is tilted to the posi¬ 
tive y direction due to the transverse field. Depending on 
the direction of the magnetization the spin Hall torque 
has a component parallel, antiparallel, or normal to the 
damping torque. This means that the work done by the 
spin Hall torque, denoted by A E s in Fig. [ 6 ] (a), is pos¬ 
itive, negative, or zero at these positions. This can be 
confirmed numerically when we calculate the work done 
by the spin Hall torque using Eq. ©• For an infinites¬ 
imal time At, the work done by the spin Hall torque 
is equal to the rate of its work ( d£ s /dt ), given in Eq. 
©, times At, i.e. A E s = (d£ s /dt)At. The solid line 
in Fig. [ 6 jb) shows an example of the calculated rate of 
the work done by the spin Hall torque (solid line), d£ s /dt 
in Eq. ©. As shown, d£ s /dt is positive, negative, and 
zero, when the magnetization undergoes one precession 
period. Similarly, the energy dissipated by the damping 
torque, d£ a /dt , can be calculated using Eq. (|5]) and is 
shown by the dotted line in Fig. [BJb). The calculated 
dissipation due to damping over a precession period is 
always negative. Details of how the rates, shown in Fig. 
[ 6 j are calculated are summarized in Appendix C. 

Note that the strength of the spin Hall torque for 
A E s > 0 is larger than that for A E s < 0 due to the an¬ 
gular dependence of the spin Hall torque, |m x (e y x m)|. 
Although it is difficult to see, thesolid line in Fig. [5(b) is 
slightly shifted upward. Thus the total energy supplied 
by the spin Hall torque during one precession, given by 
§ dt(d£ s /dt), does not average to zero and becomes posi¬ 
tive. When the current magnitude, |j|, is larger than \j c \ 
in Eq. (HOD, the energy supplied by the spin Hall torque 
overcomes the dissipation due to the damping and con¬ 
sequently the precession amplitude grows, which leads to 
the magnetization instability shown in Fig. Efa). The 
same picture is applicable when both directions of field 
and current are reversed. For this condition, the insta¬ 
bility of the magnetization is induced by the competition 
between the spin Hall torque and the damping torque. 
Therefore, the critical current depends on the damping 
constant a. When only the current direction is reversed 
in Figs. [ 6 (a) and [ 6 (b) (i.e., the sign of the magnetic field 
and current is opposite to each other), the sign of A E s is 
reversed and thus the total energy supplied by the spin 
Hall torque becomes negative. This means that the spin 
Hall torque cannot overcome the damping torque to in¬ 
duce instability. Therefore, the critical current shown in 
Eq. ([ 20 ]) only applies to the case when the sign of the 
field and current is the same. As described in Sec. m 
the same sign of the current and field in our definition 
means that the incoming electrons’ spin direction, due 
to the spin Hall effect, is opposite to the transverse field 
direction. 

Next, we consider the case when the critical current is 
damping independent. Figure El (c) schematically shows 
the precession trajectory when the applied field points to 
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(a) (b) 



FIG. 6. (a) A schematic view of the precession trajectory 

in the presence of the applied field in the positive ^-direction. 
The solid and dotted arrows indicate the directions of the 
spin Hall torque and the damping torque, respectively. The 
dashed line, which is the tangent line to the precession tra¬ 
jectory, shows the field torque. The damping torque always 
dissipates energy from the ferromagnet. On the other hand, 
the spin Hall torque supplies energy (A E s > 0) when its di¬ 
rection is anti-parallel to the damping torque, and dissipates 
energy (A E s < 0) when the direction is parallel to the damp¬ 
ing torque. When the direction of the spin Hall torque is 
orthogonal to the damping torque, the spin Hall torque does 
not change the energy (A E s s= 0). (b) Typical temporal vari¬ 
ation of the rates of the work done by the spin Hall torque, 
Eq. ©, (solid) and the dissipation due to damping, Eq. J5) 
(dotted) in the presence of the transverse field. The time is 
normalized by the period given by Eq. (IC7|) . (c), (d) Similar 
figures with the longitudinal field. 

the x direction and /3 = 0. The corresponding rate of 
work done by the spin Hall torque and the dissipation 
rate due to the damping torque are shown in Fig. [6] (d). 
Similar to the previous case, A E s can be positive, nega¬ 
tive, or zero during one precession period. However, the 
total work done by the spin Hall torque, § dt(d£ s /dt), be¬ 
comes zero in this case due to the symmetry of angular 
dependence of the spin Hall torque. This means that the 
spin Hall torque cannot compensate the damping torque, 
and thus, a steady precession assumed in the linearized 
LLG equation is not excited. This is evident in the nu¬ 
merically calculated magnetization trajectory shown in 
Fig. [5fb). For this case, the linearized LLG equation 
gives \j c \ oo, indicating that the spin Hall torque can¬ 
not destabilize the magnetization. The same picture is 
also applicable, for example, in the absence of the applied 
field and (3 = 0. 

However, an alternative mechanism can cause destabi¬ 
lization of the magnetization. As schematically shown in 
Figs. EJa) and[6l^c), there is a component of the damping 
like spin Hall torque that is orthogonal to the damping 
torque when A E s = 0. The spin Hall torque at this point 


is parallel or antiparallel to the field torque depending on 
the position of the magnetization. When the spin Hall 
torque is sufficiently larger than the field torque, the mag¬ 
netization moves from its equilibrium position even if the 
total energy supplied by the spin Hall torque is zero or 
negative. This leads to an instability that occurs before 
one precession finishes. In this case, it is expected that 
the critical current is damping-independent because the 
instability is induced as a competition between the spin 
Hall torque and the field torque, not the damping torque. 
The time evolution of the magnetization shown in Fig. 
[5] (b) represents such instability. The work reported in 
Refs . 14 i 49 discusses a similar instability condition. 

The above physical picture is also applicable in the 
presence of the fieldlike torque. The fieldlike torque, 
which acts like a torque due to the transverse field, modi¬ 
fies the equilibrium direction of the ferromagnet and thus 
the precession trajectory. Consequently, the amount of 
energy supplied by the spin Hall torque and the dissipa¬ 
tion due to damping is changed when the fieldlike torque 
is present. Depending on the sign of /?, the amount of the 
work done by the spin Hall torque increases or decreases 
compared to the case with /3 = 0. In our definition, posi¬ 
tive /3 contributes to the increase of the supplied energy, 
resulting in the reduction of the critical current. The 
complex dependence of the critical current on a arises 
when the fieldlike torque is present. 

To summarize the discussion, the critical current be¬ 
comes damping dependent when the energy supplied by 
the spin Hall torque during a precession around the equi¬ 
librium is positive. The condition that meets this criteria 
depends on the relative direction of the spin Hall torque 
and the damping torque, as briefly discussed above. To 
derive an analytical formula that describes the condition 
at which the critical current becomes damping dependent 
is not an easy task except for some limited cases^. 


VI. CONCLUSION 

In summary, we have studied the critical current 
needed to destabilize a perpendicularly magnetized fer¬ 
romagnet by the spin Hall effect. The Landau-Lifshitz- 
Gilbert (LLG) equation that includes both the damping¬ 
like and fieldlike torques associated with the spin Hall 
effect is solved numerically and analytically. The criti¬ 
cal current is found to have different dependence on the 
damping constant, i.e., the critical current scales with a 0 
(damping-independent), <a, and a 1 / 2 depending on the 
sign of the fieldlike torque. The analytical formulas of 
the damping-dependent critical current, Eqs. (HU), 
and m, are derived from the linearized LLG equation, 
which explain well the numerical results. We find that 
systems with fieldlike torque having the appropriate sign 
(/3 > 0 in our definition) are the most efficient way to re¬ 
duce the critical current. For typical material parameters 
found in experiment, the critical current can be reduced 
to the order of 10 6 A/cm 2 when ferromagnets with rea- 
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sonable parameters are used. 
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Appendix A: Initial state of the numerical 
simulation 

We assume that the magnetization in the absence of 
the applied field points to the positive z direction. In 
the presence of the field, the equilibrium direction moves 
from the z axis to the xy plane. Let us denote the zenith 
and azimuth angles of the initial state m (t = 0) as 6 and 
ip, i.e., m (t = 0) = (sin # cos (^, smOsirup, cos6). When 
the applied field points to the x-direction (n# at e x ), the 
initial state is 

{sm-yHZ wl + (pHs) 2 /H K ]\ (A1) 

V tan -1 (/3H S /H app \) J ’ 

where the value of ip is 0 < ip < 7r/2 for H a pp i > 0 and 
/3H S > 0, 7 r /2 < ip < 7 r for i7 app i < 0 and /3H S > 0, 7 r < 
ip < 37 r /2 for H app i < 0 and pH s < 0, and 37t/2 < ip < 2 tt 
for H app i > 0 and f3H s < 0. On the other hand, when 
the applied field points to the ^/-direction (n h = e y ), the 
initial state is 

[9\ f sin - 1 [(i? appl + pH s )/H K ] 

W..=H V ^ 

where the range of the inverse sine function is — 7 t /2 < 
sin -1 x < 7 t/ 2 . We note that the choice of the initial 
state does not affect the evaluation of the critical cur¬ 
rent significantly, especially in the small field and current 
regimes. 




FIG. 7. Time evolution of the z component of the mag¬ 
netization m z in the presence of the longitudinal field with 
i7 app i ~ —30 Oe, P = 0, and j m 110 x 10 6 A/cm 2 . The 
dotted line is a guide showing m z = 0. 

the critical current because the final state of the magne¬ 
tization is far away from the initial equilibrium. 

Regarding this point, we show the numerically eval¬ 
uated critical current with a different definition. The 
magnetic state can be regarded as unstable when it fi¬ 
nally arrives at a point far away from the initial stated. 
Thus, for example, one can define the critical current as 
a minimum current satisfying 

lim | m z {t) — m z ( 0)| > 5, (Bl) 

t—> oo 

where a small positive real number S is chosen to be 
0.1 here. Figure [ 8 ] summarizes the numerically evalu¬ 
ated critical current with the definition of Eq. ([Bl [) . The 
analytical formulas, Eqs. m-m, still fit well with the 
numerical results. The absolute values of the damping- 
dependent critical current are slightly changed when the 
definition of the critical current is changed. This is be¬ 
cause Eq. (1B1|) is more easily satisfied than Eq. ©, 
and thus the critical current in Fig. [ 8 ] is smaller than 
that shown in Fig. 01 However, the main results of this 
paper, such as the damping dependence of the critical 
current, are not changed by changing the definition of 
the critical current in the numerical simulations. 


Appendix B: Numerically evaluated critical current 
with different definition 

As mentioned in Sec. uni the definition of the critical 
current has arbitrariness. As an example, we show the 
time evolution of m z under the conditions of n h = e x , 
H app \ = —30 Oe, P = 0, and j = 110 x 10 6 A/cm 2 in 
Fig. [71 In this case, the magnetization initially starts at 
m z = cos[sin _ 1 (i7 ap pi/i7K)] ^0.99, and finally moves to 
a point m z 0 . 12 . Since the final state does not satisfy 
Eq. this current, j = 110 x 10 6 A/cm 2 , should be 
regarded as the current smaller than the critical current 
in Sec. m However, from the analytical point of view, 
this current can be regarded as the current larger than 


Appendix C: Energy change during a precession 

As described in Sec. El the linearized LLG equation 
assumes a steady precession of the magnetization due to 
the field torque when the current magnitude is close to 
the critical current. This is because the spin Hall torque 
compensates with the damping torque. Thus, Figs. ( 6 ^b) 
and [Hd) are obtained by substituting the solution of m 
precessing a constant energy curve of £ into Eqs. © and 
©• 

When the transverse field is applied and P = 0, i.e., 
£ = E, where E = — M f dm • H, the precession trajec¬ 
tory on the constant energy curve of E is given by^ 

rrt x (E) = (r 2 - r 3 )sn(iz, k)cn(u, fc), (Cl) 
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FIG. 8. Numerically evaluated critical currents with a different definition, Eq. dm), in the presence of (a)-(c) the longitudinal 
(uh = e x ) and (d)-(f) the transverse (n# = e y ) fields, where the value of ft is (a), (d) 0.0; (b), (e) 2.0; and (c), (f) —2.0. The 
solid lines are the analytically estimated critical current described in Sec. [Vl 


m,y(E) = r 3 + (r 2 - r 3 )sn 1 2 3 (u, k ), 


m z {E) = - r§ - {r\ - rf )sn 2 (w, k), 


(C2) T he modulus of J acobi elliptic functions is k = 
yj ( 7*2 — rs)/{ri — rs). The precession period is 


(C3) 


t{E) = 


_ 2K(fc) _ 

7v /ff appl^K/2\/ r l - r 3 ’ 


(C7) 


where u = 7 \J H { Hk/2 — r:>t , and r> are given by 


n(E) = - 


E 

M 77;ipp] 


(C4) 


r 2 (E) 


-^appl 1 

H k 



2 E 

+ MH k ' 


(C5) 


r 3 (E) 


-f^appl 



2 E 

MH k 


(C6) 


where K (k) is the first kind of complete elliptic inte¬ 
gral. The initial state is chosen to be mn y ( 0) = r%. Fig¬ 
ure [S^b) is obtained by substituting Eqs. cay dna, 
and m into Eqs. © and ©• We note that Eqs. 
(£©, and (jC3|) are functions of the energy den¬ 

sity E. Since we are interested in the instability thresh¬ 
old near the equilibrium, the value of E is chosen close 
to the minimum energy E m i n . In Fig. [6] (b), we use 
E = Emin + (E max - E m[n )/N with N = 100, where 
the minimum energy E m i n and the maximum energy 
E max are E min = -{MHk/ 2)[1 + (E app i/E K ) 2 ] and 
E max — Mif app i, respectively. The value of -f^ app i is 
chosen to be 50 Oe. Figure (6^d) is obtained in a similar 
way. 
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